Spinor condensate of Rb as a dipolar gas 
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We consider a spinor condensate of 87 Rb atoms in F = 1 hyperfme state confined in an optical 
dipole trap. Putting initially all atoms in m,F = component we find that the system evolves 
towards a state of thermal equilibrium with kinetic energy equally distributed among all magnetic 
components. We show that this process is dominated by the dipolar interaction of magnetic spins 
rather than spin mixing contact potential. Our results show that because of a dynamical separa- 
tion of magnetic components the spin mixing dynamics in S7 Rb condensate is governed by dipolar 
interaction which plays no role in a single component rubidium system in a magnetic trap. 

PACS numbers: 03.75.Mn, 05.30.Jp, 75.45. +j, 75.50.Mm 
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Since the first achievement of a Bose-Einstein conden- 
sation of trapped atomic gases, an ultracold atomic 87 Rb 
vapor is considered as the generic example of an atomic 
condensate. Rubidium atoms have many features which 
make them a perfect candidate for ultracold atoms ex- 
periments. At temperatures in the nanokelvin regime, 
interatomic interactions can be, to a very high accuracy, 
approximated by a short range two body potential with 
the s-wave scattering length a s being the only relevant 
parameter. In a magnetic trap, the two low field seek- 
ing states stable with respect to two body processes, are 
the two hyperfme states of the atomic rubidium ground 
state: \F = 2,mp = 2) or \F — \,nip — —1). Simple 
estimation of a characteristic contact energy gives: 



E c = (ATih 2 a s /M) n . 



(1) 



where a s is the s-wave scattering length, M is the mass 
of the atom and n is a typical atomic density. The spin 
dipole-dipole interaction energy is of the order of: 



E d = fi n, 



(2) 



where magnetic moment of 87 Rb in an F = 1 hyperfme 
state is equal to (i = \[Ib where \ib is the Bohr magne- 
ton. The ratio of these two energies is: 



E d /E c = 4.2 x 10- 



(3) 



Due to the smallness of the atomic magnetic moment the 
dipole-dipole interaction between atoms is several orders 
of magnitude smaller than the contact interaction and 
can be ignored. The rubidium condensate in a magnetic 
trap is thus a perfect example of a weakly interacting 
system with contact interaction characterized by a single 
parameter - the scattering length. 

Experimental achievement of a condensate of atoms 
with sizeable long range dipolar forces was a great chal- 
lenge. Finally, after overcoming a number of serious dif- 
ficulties a chromium 52 Cr condensate characterized by 
a relatively large dipole moment was created More- 
over, by utilizing the technique of Feshbach resonance the 



contact term was practically turned off and pure dipolar 
condensate was achieved [2|. 

In this paper we show that yet another and experimen- 
tally much simpler way to dipolar condensates is possi- 
ble. Namely, optical dipole traps allow for a simultaneous 
trapping of various magnetic components of a given hy- 
perfme state, i.e. spinor condensates. A spin dynamics 
of rubidium F = 1 and F — 2 states and the formation 
of a condensate due to increase of the atomic number in 
a given Zeeman sublevel was studied experimentally 0|. 
In the experiment Q a transfer of atoms from the initial 
roj? = state to tuf = ±1 states was observed. 

The theoretical studies related to \4i were performed in 
one or two spatial dimensions 0, 0,13| and dipole-dipole 
interactions so far were ignored. On the other hand, some 
theoretical and experimental evidence of enhanced role 
of dipolar forces in F = 1 87 Rb was reported. This is 
the observation of disintegration of a helical structure of 
magnetization d| or Einstein-de Haas effect Some 
authors @, [HI, E4] have already stressed an enhanced role 
of dipole-dipole interactions in a spin dynamics of ferro- 
magnetic 87 Rb. This is because the contact spin mixing 
term (see the Hamiltonian Eq. ©) is proportional to 
E a = Aith 2 {a,2 — ao)n/3M where the scattering lengths 
ao = 5.387nm and £12 = 5.313nm determine collisional 
cross sections in a channel of a total spin and 2 respec- 
tively. Therefore a ratio of the dipolar energy to the spin 
mixing contact term in the ferromagnetic rubidium is as 
large as: 



E d /E s = 0.09. 



(4) 



In order to analyze processes responsible for a spin dy- 
namics we shall discuss a role of different terms of the 
Hamiltonian of the system: 



H = d 3 r (£t( r )ff £.( r ) - 7 *t (r)BF ..j.( r) 



-c *t( r )|r|(r)^(r)^.(r) 
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d 3 rd 3 /*t(r)$t( r ')^ w(r _ r ')^.( r ')$^( r ) 

(5) 



where repeated indices (taking values +1,0 and -1) are 
to be summed over. 

The first term in J5]) is the single particle kinetic 
energy and the trapping potential energy V tr : H = 
-|gV 2 + Vt r (r). The second term describes the inter- 
action with magnetic field B with 7 being the gyromag- 
netic coefficient which relates the effective magnetic mo- 
ment with the hyperfine angular momentum (fj, = 7F). 
The terms with coefficients cq = A-Kh 2 {aQ + la^jZM and 
C2 = 4irh 2 (a2 — ao)/3M describe the spin-independent 
and spin-dependent parts of the contact interactions. Be- 
cause C2 < the ground state of 87 Rb is ferromagnetic. 
F are the spin-1 matrices. The last term describes the 
magnetic dipolar interaction of two magnetic dipole mo- 
ments located at r and r': 
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|r-r'| |r-r'|" 
(BV(r-r'))(F H .(r-r')) . (6) 



The field operator ^(r) annihilates an atom in the 
hyperfine state \F = 1 , i) at a point r. Using the classical 
fields approximation [13j we replace the field operators 
^(r) by the classical wavefunctions ^j(r). The equation 
of motion for these wavefunctions is 



ih 



d_ 

"St 



*1 



= {H + H B +H C + H A ) 



•(7) 



The diagonal part of H c is given by H c \\ — (co + 
ca)***i + (co + c 2 )*5*o + (co - ca)**!*-!, # c oo = 

(C + C 2 )*J* X + C *o*0 + ( c + C 2 )** flc-1-1 = 

(c - + (c + c 2 )^ + (c + cajtfii*.!. The 

off-diagonal elements that describe collisions not pre- 
serving the spin projection of each atom are equal to 
H cW = C2**i* 0) H cO-i =c 2 ^V Moreover H c i_i = 
0. On the other hand, for the 2?^ term one has H^j = 
/ dV**(r')T^ (^(r'). This term is responsible for the 
change of total spin projection of colliding atoms. 

In the following we use the oscillator units where a 

1 II 

distance is measured in aho — (h/Mui) where u = 
2ir x 100 Hz. Initially the system is excited with all 
atoms in mp = component (with mp — ±1 compo- 
nents equal to zero). As the first step we compute the 
mp = ground state of the system by means of the imag- 
inary time propagation. The resulting wave function was 
then randomly perturbed in order to inject about 10% of 
the excess energy. 

The details of the classical field approximation are re- 
viewed in [13]. Within the classical field approxima- 
tion this amounts to introduce a thermal fraction into 
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FIG. 1: (color online) Populations of thermal clouds for tuf = 
+1 (red), niF = (black), and mf = -1 (green) states as a 
function of time with (the main frame) and without (the inset) 
dipole-dipole interactions. The parameters are N=3xl0 5 , /3 = 
1. Total populations of trf = ±1 components are identical if 
dipole-dipole interaction is turned off. The external magnetic 
field is equal to zero. 
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FIG. 2: (color online) Kinetic energy oimp = +1 (red), tuf = 
(black), tuf = — 1 (green) states as a function of time with 
dipole interactions (upper panel) and without them (lower 
panel) for N = 3x10° and /3 — 1. The time of thermalization 
is shorter if dipolar interactions are included. The external 
magnetic field is equal to zero. 



the initial state. In all our simulations 20% of atoms 
are in the thermal clouds, corresponding to T/T c w 0.58. 
Moreover, we put a small seed (0.3% of all atoms) into 
mp = ±1 components what is necessary to initiate a spin 
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dynamics. We study numerically the dynamics leading to 
a thermal equilibrium. In all simulations we use a grid 
of 42 points in each direction with spatial step equal to 
Sx = 0.6. 

We start our study with a spherically symmetric sys- 
tem, j3 = u> z /uj r = 1 (where u> r and uj z are the radial 
and axial trap frequencies respectively) of N = 3 x 10 5 
atoms, uj r = 2tt x 100 Hz and B = 0. Due to interac- 
tions the initially empty magnetic components become 
populated and finally all three thermal clouds oscillate 
around the same value what signifies a thermal equilib- 
rium, Fig. [TJ At the equilibrium, populations of ther- 
mal fractions of mf = ±1 components fluctuate inde- 
pendently. Moreover the kinetic energies accumulated 
in every component saturate and equalize, Fig. [21 The 
total populations of mf — ±1 components need not be 
identical in the presence of dipolar interactions due to 
a noise triggered spontaneous chiral symmetry breaking. 
To our surprise the dipole-dipole interactions, seemingly 
about order of magnitude smaller than contact spin mix- 
ing term, significantly decrease the thermalization time: 
from 1.2s without dipole-dipole term to 0.35s when this 
long range term is retained. 

In order to get a better insight into the origin of the 
observed enhanced role of the dipole-dipole interactions 
we examine in more detail a spatial structure of differ- 
ent magnetic components. Typical density profiles are 
shown in Fig. [3j At t = 250ms (upper panel) densities 
and phases (not shown) of multicomponent spinor wave 
function indicate the existence of a coreless vortex with 
winding numbers equal (—1, 0, +1) for tuf = (+1, 0, — 1) 
respectively. This vortex disappears on the time scale 
of milliseconds. Let us note that density profiles in 
mp = +1 and mp = —1 components are not axially 
symmetric and are rotated with respect to each other 
by 7r/2. The different magnetic components have small 
spatial overlap. Similar separation of magnetic phases 
persists during the evolution, see Fig. [3] (lower panel). 
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FIG. 3: The atomic density at z = plane of tuf = +1 
(left frame), tuf = (middle frame) and nif = -1 (right 
frame) components for the spherically symmetric trap and 
N = 3 x 10 5 atoms at time t = 0.25s (top) and t = 0.95s 
(bottom). The external magnetic field is equal to zero. 
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FIG. 4: Ratio of averaged dipole-dipole energy to the mean 
value of the contact spin mixing term £d/e c as a function of 
time. The external magnetic field is equal to zero. 
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FIG. 5: (color online) Kinetic energy of mp = +1, 0, —1 com- 
ponents for TV = 10 5 atoms as a function of time. Top figure 
- the oblate geometry, (3 = 2, while the lower one - the pro- 
late shape, (3 = 0.5. Results were obtained from dynamics 
generated by the full Hamiltonian - main frame, and by the 
Hamiltonian without the dipole-dipole term - insets. The ex- 
ternal magnetic field is equal to zero. 

The transfer of atoms to initially empty magnetic com- 
ponents arises through dynamical instability. At ini- 
tial stages of evolution the contact spin mixing term 
is responsible for the dynamical instability leading to 
the formation of domains of opposite magnetization 
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[ana iia .[l7|. which only weakly overlap. The dynam- 
ically formed spinor wave function cannot be described 
by a single mode approximation models [ill. fl2l. [l8| . The 
simplest estimation of a magnitude of the contact spin 
mixing term Eq. U) and the dipole-dipole energy Eq. 
makes use of a mean atomic density. This estimation is 
not adequate if the magnetic domains are formed because 
the dipole-dipole interaction is the long range one as op- 
posed to the zero range contact potential. When spin 
domains are formed the role of the contact term signif- 
icantly decreases and dipole-dipole interactions become 
dominant if the spin dynamics is concerned. A mean 
value of the contact spin-mixing term averaged over the 
spatial distribution of the spinor wave function is: 

e c = c 2 J d 3 rn(rW^(r)^ (^o(r), (8) 

while the dipole-dipole averaged energy reads: 

e d = -h^ 2 J d\ J rf 3 ^(r)^$ (r) , (9) 

where V{v,v') = 3/\/2e-^cosesine(|1' 1 (y)| 2 - 
|*_i(r')| 2 ) + 3/2e- 2 ^sin 2 e(*t(r')* (r') + 

#S(r')*-i(r')) (1 3/2sin 2 e)(^(r')*i(i-') + 
*I 1 (r')\E , (r')), and 0, 9 are the spherical angles of 
R = r — r'. The ratio of these two terms is shown in 
Fig. 21 During the evolution the dipole-dipole term 
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FIG. 6: Densities of nif = +1 (left frame), niF = (middle 
frame) and tuf = — 1 (right frame) components for the case 
of N = 10 5 atoms. The top frame is for /3 = 2 at t = 0.7s, and 
the bottom frame is for (3 = 0.5 at t = 0.35s. The densities of 
my = +1 and m_F = — 1 components for (3 = 2 are spatially 
separated while for /3 = 0.5 they overlap. 

is comparable to the contact one and frequently the 
dipolar interaction becomes dominant. This is the main 
observation of our study. 

The structure of the dynamically unstable modes de- 
pends on a geometry of the system. To explore this fact 
we analyzed two cases, j3 = 2 (u> z — 2ir x 200 Hz, u> r = 2tt 
x 100 Hz) corresponding to an oblate profile and (3 = 0.5 
(tu z = 2ir x 100 Hz and ui r — 2ir x 200 Hz) characteristic 



for a prolate shape. In the case of the oblate geometry 
the dipolar forces significantly accelerate the spin dynam- 
ics while no effect is observed for the prolate geometry, 
Fig. [5j The spatial density profiles in Fig. [6] help to 
understand this difference. For the oblate shape the dy- 
namics leads to the separation of phases and formation 
of only weakly overlapping domains of opposite magne- 
tization. This is the reason why the role of zero range 
interaction terms significantly decreases and long range 
dipole forces start to dominate the spin dynamics. The 
same dynamics in the case of the prolate geometry favors 
the formation of almost identical structures of = 1 
and rriF = — 1 components and spin mixing contact term 
dominates over the dipole-dipole interaction. 

The theoretical studies of dipole-dipole interactions in 
F = 1 87 Rb [H] show, that it is not easy to observe dipo- 
lar effects at non-zero magnetic field. In particular, dipo- 
lar effects are easy to observe at magnetic field ~ 10/iG 
[lfj| . Our study shows, that it is experimentally possible 
to see dipolar effects by observing time of thermalization 
of the system. At T — OK, the time of thermalization 
for both type of systems (with and without dipole-dipole 
interactions) is the same and equals t w 0.75s. For non- 
zero temperatures (Nq/N ~ 0.8, where Nq is a number of 
condensed atoms and N is the total number of atoms) the 
time of thermalization is different for both cases. The 
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FIG. 7: (color online) Populations of thermal clouds for niF = 
+1 (red), uif = (black), and tuf = — 1 (green) states as a 
function of time with (the main frame) and without (the inset) 
dipole-dipole interactions. The parameters are N=3xl0 5 , /3 = 
1, and B z — lmG. The time of thermalization is t ~ 0.4s. 



time of thermal equilibrium at B z = lmG (Fig. [7]) for 
spherical trap is almost identical as for B z = (Fig. [TJ. 
The time of thermalization is equal t w 0.4s and is still 
three times smaller with the dipole-dipole interactions 
included. For B z = lOOmG, the system is thermalized 
after t ~ 0.7s (Fig. [TT]). For the without dipolar in- 
teractions the system is thermalized after t » 1.2s (this 
time is the same as without external magnetic field). So, 
even at higher magnetic fields, we can see dipolar effects 
by observing the time of thermalization of the system. 
The spatial density profiles (Fig. [9]) show, that in this 
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FIG. 8: (color online) Kinetic energy of mp = +1 (red), 
rriF = (black), m,F — —1 (green) states as a function of 
time with dipole interactions (the main frame) and without 
them (the inset) for spherical symmetry. The time of thermal- 
ization at B z = lmG remains shorter if dipolar interactions 
are included. 
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FIG. 9: (color online) The atomic density at z = plane of 
m,F = +1 (left frame), itif = (middle frame) and tuf = — 1 
(right frame) components for the spherically symmetric trap 
and N = 3 x 10 5 atoms at B z — lmG at time t = 0.45s with 
(top) and without (bottom) dipole-dipole interactions. 



FIG. 10: (color online) The sequence of frames representing 
changing of spin textures with dipole-dipole interactions. The 
value of external magnetic field is B z = lmG and Larmor 
precession frequency is <^i'" G = 6.9 and it = 1.45 x 10' 3 s. 
The sequence of frames starts from 1.2s and ends after the 
single cycle. 
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relatively strong magnetic field the unstable Bogoliubov 
modes have m = ±1 azimuthal quantum number in a 
close analogy with the 2D simulations of The latter 
had no dipole-dipole forces, thus we conclude that the 
magnetic field tends to mask the dipolar effects [l9| . 

The dynamics of both systems allows to observe rich 
structures of spin textures. FigurefTUlshows typical struc- 
tures of [F x ,F y ] quantity for dipolar case at B z = lmG 
(Fig. OH and B z = lOOmG (Fig. [TJJ|. 

For non-dipolar case we can observe similar structures 
in 'xy' plane. Additionally, rotation of spin domains is 
observed. During a single rotation, it is possible to ob- 
serve different structures of spin textures. Using Larmor 
formula 

u L = lB z , (10) 



FIG. 11: (color online) Populations of thermal clouds for 
trif = +1 (red), tuf = (black), and tuf = — 1 (green) states 
as a function of time with (the main frame) and without (the 
inset) dipole-dipole interactions for spherical symmetry and 
for B z — lOOmG. The time of thermalization at B z = lOOmG 
is longer then for B z = lmG, but still remains shorter if dipo- 
lar interactions are included. The system thermalizes after 
t w 0.7s. 



where u>l is the Larmor frequency we can calculate 
the number of rotations per unit time. For example, 
w imG = 6 9 and t v,nG = 1.45 x iq-3 s . For B z = lOOmG, 

w 100mG = g 90 and t 100mG = 145 x 10 -5 g which is con . 

firmed by our simulations. For low magnetic field we 
observe a different time of spin precession. For exam- 
ple, at B z = 10/iG the rotation period for dipole-dipole 
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FIG. 12: (color online) Kinetic energy of tuf = +1 (red), 
m,F = (black), m,F — —1 (green) states as a function of 
time with dipole interactions (the main frame) and without 
them (the inset) for spherical symmetry for magnetic field 
B z — lOOmG. The system thermalizes after t ^ 0.7s. 
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FIG. 13: (color online) The sequence of frames representing 
changing of spin textures. The value of external magnetic field 
is B z = lOOmG and Larmor precession is u^ 00mG = 690 and 
ii = 1.45 x 10~ 5 s. The sequence of frames starts from 0.75s 
and ends after the single cycle. The dipole-dipole interactions 
are included. 



interaction case is about two times smaller. It can be 
explained by replacement of external magnetic field by 
effective magnetic field which includes external magnetic 
field, the mean field originating from the spin-exchange 
and the mean field originating from the dipole-dipole in- 
teractions [20] . 

In summary we have shown, that the dipole-dipole in- 
teractions play a major role in the dynamics of the 87 Rb 
condensate in F = 1 hyperfine state: a) the spin chang- 
ing s-wave scattering length is small (noted before), b) 



the contact interactions trigger the formation of weakly 
overlapping magnetic domains, c) because of a) and b) 
the importance of contact interactions with respect to 
the long range forces is reduced, d) as a consequence 
the thermalization time is significantly shorter with than 
without the dipolar term. We hope, that this findings 
may be verified experimentally by monitoring dynamics 
of the thermal clouds of all magnetic components turning 
off and on the dipole-dipole interactions 
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